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Recently, unusual and strikingly beautiful seahorse-like 
growth patterns have been observed under conditions of quasi- 
two-dimensional growth. These 'S'-shaped patterns strongly 
break two-dimensional inversion symmetry; however such bro- 
ken symmetry occurs only at the level of overall morphology, 
as the clusters are formed from achiral molecules with an achi- 
ral unit cell. Here we describe a mechanism which gives rise to 
chiral growth morphologies without invoking microscopic chi- 
rality. This mechanism involves trapped electrostatic charge 
on the growing cluster, and the enhancement of growth in re- 
gions of large electric field. We illustrate the mechanism with 
a tree growth model, with a continuum model for the motion 
of the one- dimensional boundary, and with microscopic Monte 
Carlo simulations. Our most dramatic results are found using 
the continuum model, which strongly exhibits spontaneous 
chiral symmetry breaking, and in particular finned 'S' shapes 
like those seen in the experiments. 

PACS numbers: 47.54.+r,61.43.Hv,68.70.+w 



I. INTRODUCTION 

Growth phenomena are known for the complexity and 
beauty of the patterns they can lead to (lj . Most of this 
complexity results from different kinds of instabilities as- 
sociated with growth, such as the Mullins-Sekerka insta- 
bility of growth fronts ||, or fingering instability ||-||. 
The presence of instabilities implies that a tiny micro- 
scopic noise can result in macroscopic changes of shape, 
and hence lead to a variety of shapes. (For instance, the 
formation of snowflakes, a growth phenomenon familiar 
to everyone, produces hundreds of different shapes ||). 
However, despite the variety of shapes, the vast major- 
ity of growth patterns preserve left-right symmetry; in 
other words, essentially all of the known growth patterns 
are achiral. One well-known and historic exception is 
the appearance of hemihcdral faces on crystals, yielding 
facetted forms which are not invariant under inversion 
0. In the mid-19th century, such faces were identified 
in molecular crystals of sodium ammonium tartrate by 
Pasteur, and the broken chiral symmetry was ascribed 
by him to the microscopic chirality of the constituent 
molecules. 



In this work we concentrate on two-dimensional 
(or quasi- two-dimensional) growth forms. For two- 
dimensional forms the relevant inversion operator is x — > 
—x or y — > —y, but not both; we will call such an opera- 
tion "2D inversion" , and forms distinguishable from their 
2D inverse "2D chiral". (Also, since we concentrate en- 
tirely on 2D henceforth, we will sometimes shorten these 
terms by omitting the leading "2D" qualifier.) 2D chiral 
growth forms arc not common In those rare cases 
where chiral growth patterns do appear &|l(| , the inver- 
sion symmetry is already broken at the microscopic level. 
An example is the formation of spiral crystals B during 
the compression of a phospholipid monolayer on a water- 
air interface. In this case the individual phospholipid 
molecules possess a 3D chirality. Each molecule also has 
a preferential orientation (hydrophilic head down) with 
respect to the water-air interface. This consistent orien- 
tation then gives a monolayer with two-dimensional in- 
version symmetry already broken |p"l| at the microscopic 
level — assuming only that the molecules themselves have 
a predominance of a single enantiomer. And in fact, for 
spiral crystals to appear, one needs to have a monolayer 
consisting predominantly of a single enantiomer. The 
handedness of the crystals depends directly on the hand- 
edness of the dominant enantiomer, and no chirality ap- 
pears for racemic monolayers || . While there are several 
competing explanations of how the microscopic chiral- 
ity leads to the macroscopic chirality Jl^,|l^], it is nev- 
ertheless clear that the latter occurs only because of the 
former. Similarly, in another chiral growth example — the 
formation of chiral bacterial colonies jl(J — the individual 
particles (bacteria) also have a three-dimensional chiral- 
ity (of a single 'sign') which then manifests itself as a 2D 
chirality when coupled with a 2D substrate. The bacte- 
rial aggregates are observed to be (2D) chiral, and always 
with the same handedness jtO|. Thus, in each of these 
cases, it is clear that the macroscopic 2D chirality of the 
aggregates results from a microscopic 3D chirality of the 
elementary building blocks. 

Recently, a novel and very beautiful type of growth 
pattern has been reported jljjl^]. A typical pattern 
strikingly resembles a seahorse (in the form of an 'S' 
shape, with 'fins' on the outer curved edges), and so 
has a strongly broken 2D inversion symmetry. The pat- 
terns were discovered during growth studies of fullerene- 
tetracyanoquinodimethane (Ceo-TCNQ) thin films. Sub- 



1 



sequently, very similar patterns were obtained using 
TCNQ only fig] . The broken symmetry is one of the 
most striking aspects of the patterns, as well as one of the 
principal mysteries connected with them. The mystery 
arises because — in contrast to the two cases mentioned 
above — in these experiments there is no microscopic sym- 
metry breaking: TCNQ molecules are themselves inver- 
sion symmetric [fl7f . Furthermore, even though the 'sea- 
horse' aggregates are polycrystalline Q, one can prob- 
ably rule out symmetry breaking at the level of the unit 
cell, since TCNQ crystals are also achiral p7| , p!8[ . 

It is however important to note that both left- and 
right- handed patterns appear in approximately equal 
numbers Jlq] , Thus, on average, the experiment does 
not break inversion symmetry; instead the symmetry is 
broken spontaneously, for each island, during the growth. 
That is, the "seahorse" growth experiments represent an 
almost unique case of spontaneous 2D chiral symmetry 
breaking during quasi- two-dimensional growth. 

We say 'almost' unique because we are aware of only 
one other growth phenomenon exhibiting such sponta- 
neous symmetry breaking, namely, phyllotaxis: the pat- 
tern of leaves, buds, scales, etc in growing plants pgj . 
As demonstrated by Douady and Couder this phe- 
nomenon can be understood cleanly in two dimensions; 
and furthermore, the resulting spiral growth patterns are 
clearly chiral, and the symmetry breaking is clearly spon- 
taneous ■ Outside of this one example from botany, 
however, we know of no example of two-dimensional 
growth — experimental or theoretical — in which the re- 
sulting growth patterns spontaneously break 2D inver- 
sion symmetry. 

In this work we construct a growth model which does 
yield such spontaneous symmetry breaking. More pre- 
cisely, we consider a set of models, all embodying the 
same ideas. These ideas involve a novel form of long- 
ranged branch competition and growth, arising from elec- 
trostatic effects. We have found that such a mechanism 
can lead to growth forms which spontaneously break two- 
dimensional inversion symmetry. 

The growth models which we will consider share several 
important properties. These properties are simple, and 
can be formulated independently of the nature of the un- 
derlying physical processes. The physical picture which 
we consider involves the following elements: (i) branch- 
ing, that is, every growing branch should eventually give 
rise to new branches; (ii) strong branch competition — in 
fact, the competition has to be so strong, that only 2 
main branches "survive" ; (iii) long range branch repul- 
sion: the two branches need to "feel" one another and 
curve away from each other. 

Branching is a very common property in growth phe- 
nomena Q. Branch competition is also very common. 
It is usually caused by screening — that is, by the com- 
petition between growing branches for incoming parti- 
cles. But competition due to screening alone is not strong 
enough to lead to two- armed shapes. For instance, the 
diffusion-limited aggregation (DLA) model leads to clus- 



ters having 4 or more branches |21|-|25[]. Competition 
for incoming particles may also cause some branch re- 
pulsion. But this effect is obviously very short-ranged 
(a branch "feels" only its neighbors). Thus in order to 
achieve (i)-(iii), one needs to introduce a long-ranged in- 
teraction into the system. As will be discussed later in 
section III A , electrostatic forces may play an important 



role in the formation of the seahorse patterns. Hence, in 
our models, the long range interaction between branches 
is also of electrostatic origin. 

We will start with a very simple deterministic "tree" 
growth model which has properties (i)-(iii) by construc- 
tion, and show that this model prefers chiral rather than 
symmetric shapes. We will then consider a more realistic 
quasi-equilibrium continuum model, in which properties 
(ii)-(iii) arise naturally due to electrostatic interactions. 
This model yields two-armed, finned, S-shaped forms for 
a range of growth parameters, for essentially the same 
physical reasons as does the tree model. We also re- 
port some preliminary studies involving the same physi- 
cal ideas but using a microscopic Monte-Carlo approach. 
These modeling efforts are inspired by the puzzling and 
remarkable experimental patterns; and they yield quali- 
tatively similar growth forms. It is also encouraging that 
further growth experiments involving a static in-plane 
electric field (which we discuss briefly below, and in de- 
tail in another paper jl6|]) have provided support for our 
ideas. 



II. CHARGED TREES 

In this section we demonstrate, using a simple and 
highly schematic model, how long-ranged branch repul- 
sion and competition may cause chiral symmetry break- 
ing. We will formulate a simple growth model where such 
branch repulsion and competition are present by con- 
struction, and show that chiral 'S' shapes are preferred 
energetically over symmetric shapes. 

A charged tree model is constructed as follows (Fig. [l]). 
A tree starts as as a single charged rod. The ends of this 
rod are considered to be "alive" . Then each alive branch 
emits two branches: one branch of length and the other 
of length l\. Both new branches grow at a predefined 
angle 0. All three quantities — Iq, h, and 9 — are the same 
for opposite ends, and do not vary during the growth. 
When lj / Ii, there are 4 possible combinations of growth 
on every step; among these, the model chooses the tree 
with the lowest electrostatic energy. Then the longer of 
the newly added branches become new "alive" branches, 
the shorter ones "die", and the process is repeated again. 
If two (or more) configurations have the same energy then 
the selection is done randomly. The first 2 steps of growth 
are shown in Fig. |l|. 

To complete this model we need to specify how we will 
compute the electrostatic energy, as the energy of a ID 
charged rod diverges. The most obvious way to deal with 
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this problem is to assign some small (but finite) width w 
to the branches, with this width satisfying w << lo,h- 
This allows us to work with 2D, rather than ID, charge 
density. To compute the electrostatic energy of the tree, 
we further break the branches into smaller (linear) pieces. 
The ith piece has the length i,, and the linear charge den- 
sity Xi, which is in an obvious way related to 2D charge 
density = Xi/w. 

We will consider 2 possible charge distributions: f) a 
conducting charge distribution f7, = const, where Ui is 
the potential at the center of the i-th piece; 2) a uniform 
charge distribution <jj = const. In both cases the total 
charge of the system will be normalized by the require- 
ment that the average linear charge density be equal to 
unity, that is 



A 



J2 = 1 



(2.1) 



where L = ^2U is the total length of the structure. 

i 

To find the charge densities for each piece we construct 
a set of linear algebraic equations 
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(2.2) 



where 



i = 0, . . . , N — 1, and f/y is the electrostatic po- 
tential which would be induced at the center of the z-th 
piece by the j-th piece if the latter had a u nit charge 
density. Clearly the solution {A^} of Eq. fl2.2| ) meets the 
requirement Ui = const (the constant was set to 1); this 
solution can then be e asily rescaled to meet the normal- 
ization condition (2.1). 

Using elementary electrostatics, one can show that the 
constants Uij are given by 



,N 



U u = 2m(5.44tu/y 



(2.3) 



and for i ^ j 



^/(l j /2 + x i r + yf-(l j /2 + x i ) 

where lj is the length of j-th piece, Xi , yi are the coordi- 
nates of the center of the «-th piece (the origin is assumed 
to be at the center of the j-th piece, and the x-axis di- 
rected along the j-th piece). 

The growth of such trees was studied numerically, and 
all trees grown according to these rules demonstrate some 
chirality. One of the typical S-like shapes is shown in 
Fig. | 

Why does the tree prefer to break the left/right sym- 
metry? To answer this question, we note that, by intro- 
ducing the repulsive electrostatic interaction, we effec- 
tively made the tree keep its branches as far away from 
each other as possible. This observation alone accounts 
for the fact that the tree "chose" the 'S' shape on the 



first split. However, during subsequent splits the result 
is determined by the interaction between all branches, 
and the outcome depends crucially on how the charge 
is distributed over the tree: the more charge is concen- 
trated on the small "dead" branches, the stronger is the 
symmetry breaking (i.e. the higher is the achievable cur- 
vature of the main arms). We can illustrate this point 
by considering the case of a uniform charge distribution 
Xi = const. In this case the potential energy of a tree is 
given by 



1 A 2 

^ u = 9 S UijXjXik = ) UijU 



(2.5) 



We have grown a number of trees using the above rules, 
but with the assumption of a uniform charge distribution. 
The resulting trees have a much more weakly broken sym- 
metry, which comes primarily from the symmetry break- 
ing at the first branching; and they are not S-shaped. 

This result is consistent with the idea that a higher 
charge density on the external "dead" branches (which 
results for instance from the conducting charge distribu- 
tion) enhances the overall chirality of the cluster. We 
will return to the question of how the charge redistribu- 
tion influences which branches survive or die in the next 
section. 

We also note that one can construct a non- 
deterministic charged tree growth model, where the 
growth rates themselves are determined by the elec- 
trostatic interaction between branches. This non- 
deterministic growth model also leads to a spontaneous 
chiral symmetry breaking p6[ . 



III. THE CONTINUUM MODEL 

A. Construction of the Model 

In this section we will consider a more realistic growth 
model, in which strong branch competition and survival 
of only 2 main arms occur naturally (for some range of 
parameters). That is, for the model we now describe, 
properties (ii)-(iii) (strong branch competition and re- 
pulsion) arise as a result of the dynamics of growth. In 
common with the tree models described above, a crucial 
ingredient is the presence of a long-ranged electrostatic 
interaction. 

First we will consider in more detail what actually hap- 
pens in the "seahorse" experiments [ p^i|JT5[ . In these ex- 
periments layers of TCNQ are deposited using the ion- 
ized cluster beam (ICB) deposition [^,^8| method. With 
this method the TCNQ molecules are ionized and then 
accelerated towards the substrate, where they arrive with 
high kinetic energy (and thus high mobility) and there- 
fore can diffuse along the substrate and form growing 
clusters. A small fraction a of the diffusing particles 
(estimated jl6| to be typically ~ 10%) are charged (all 
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of the same sign). Thus the growing islands will also 
carry some (time-dependent) charge. The magnitude and 
time dependence of this charge are not known; it depends 
on many complicating factors, including the repulsion of 
the charged, diffusing particles by the charge on the is- 
land; leakage to the substrate, both from the diffusing 
charged walkers and from the charged aggregate; and 
the "rain" of charged particles directly on the growing 
islands. We will treat this time-dependent charge Q(t) 
in an extremely simple way below; our motivation is to 
explore the kinds of effects that electrostatic charge may 
have on growth processes. 

The field of the island, in contrast with the field of 
walkers, is not random, and therefore will play the dom- 
inant role in how the particles (walkers) diffuse. Hence 
in our work we neglect the random field of the walkers. 
We assume that neutral walkers have a non-zero polar- 
izability; hence the diffusive motion of both charged and 
neutral walkers is affected by the electrostatic field of the 
aggregates. To consider the simplest case, we will neglect 
any possible cluster-cluster interactions and consider only 
the growth of an isolated island. 

An obvious approach to this problem would be a Monte 
Carlo (MC) simulation. However, a simple estimate of 
the number of particles in seahorse clusters in the ex- 
periments gives f0 6 -10 7 particles. A direct microscopic 
Monte-Carlo simulation for a problem of this size is very 
hard, if not impossible. We have performed some limited 
MC studies (involving much smaller particle numbers, ie 
N ~ fO 3 ), which we will describe briefly in section [TV| . 
Here we will consider an alternative approach, in which, 
instead of tracing the motion of individual particles, we 
will compute local growth rates for the cluster boundary, 
which we treat as a continuous ID curve. We will obtain 
equations for the motion of this f D curve, and study the 
kinds of growth that result. 

First, let us consider a growing island surrounded 
by diffusing walkers. If a walker hits the island, then 
with some probability p s ( "sticking" probability) it (the 
walker) becomes a part of the island. Then the local 
growth rate is 



G(x,y,t) 



dh(x,y) 
dt 



h(x, y)5; 



d 2 N hlts {x,y) 
dtdl 



Ps (3.1) 



where dh(x, y) is the displacement of the given boundary 
point (x, y) during the time dt, h(x, y) is the unit vector 
normal to the boundary, S m is a typical intermolecular 
distance in the growing cluster (thus, is the area oc- 
cupied by a single molecule), and d 2 Nhits{x,y)/{dtdl) is 
the number of hits per unit boundary length per unit 
time. We then take 



d 2 N hits (x,y) 
dtdl 



= a hlts N(x,y)v T 



(3.2) 



where N(x, y) is the walker concentration near the point 
(x, y) on the boundary, vt is an average thermal velocity 
of the walkers, and ahits is a numerical factor of order of 



unity. Thus, we can rewrite the equation for local growth 
rates as 



G(x,y,t) 



dh(x,y) 
dt 



n{x,y)G T p s N(x,y) (3.3) 



where Gt = ct/iits^T^m is a constant which depends only 
on the temperature. 

Here we will assume that the growth is slow enough 
to be considered as a quasi-equilibrium process (which 
is the case when p s <c 1). Our motivations for this as- 
sumption are twofold: first, it is physically motivated, 
in that the sticking probability may indeed by very 
small, due to the high kinetic energy of the walkers in 
an ICB experiment; and second, it makes the problem 
tractable, giving a simple form for N(x, y) which en- 
ables us to concentrate on the motion of the boundary. 
Given this assumption, then, the concentration of walk- 
ers is given by a quasi-equilibrium Boltzmann distribu- 
tion: N(x,y,t) = M ) exp[~U(x,y,t)/kT] where U(x,y) 
is the potential energy of a walker at the point (x, y), 
is the concentration far away from the cluster, and k is 
Boltzmann's constant. 

In our case there are 2 different kinds of particles 
present in the system: charged and neutral walkers, each 
kind having a different concentration and a different po- 
tential energy function. The overall growth rate is given 
by 



dh(x, y, t) j dt = dh n /dt + dh c /dt 



(3.4) 



where the subscript n denotes neutral walkers, and the 
subscript c charged walkers. The walker concentrations 
are given by 

Np{x,y,t) = Nf ] exp[-Up(x,y,t)/kT] (3.5) 



with (3 = n or c and (as discussed above) Nc <C Nn 
The potential energy U(x, y) is equal to —\E 2 /2 for neu- 
tral walkers (of polarizability x) hi a field E = \E(x,y,t)\, 
and V(x, y,t)e for charged walkers. 

Now we assume that the cluster is conducting. As 
with the quasiequilibrium assumption, our reasons are 
both computational and physical: the conducting cluster 
is rather straightforward to treat numerically (and even 
analytically in some special cases); but also, from our 
charged-tree studies, we expect that a conducting cluster 
will enhance the type of branch competition which we 
wish to study here. For a conducting cluster, U c (and 
hence dh c /dt) are each independent of position on the 
island. The applicability of this assumption to the sea- 
horse experiments will be discussed below. 

The model we have constructed thus far requires the 
computation of the electric field due to the charge on a 
two-dimensional growing cluster, which is in general of 
an irregular shape. This electric field is determined by 
the charge distribution <j(x, y, i) on the island. However, 
E(x, y) diverges near the edge of any 2D charge distri- 
bution; hence, instead of using the field at the edge, we 
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will use the field at a small (molecular) distance from the 
edge. That is, instead of E(x, y, t) = E(f, t) we will use 
E(r + 8 ■ h(x,y),t), where 8 is the 'sticking distance' — 
the distance at which a diffusing particle sticks to the 
cluster; we will assume 8 ~ 8 m . The electric field is de- 
termined by the shape of the island, to within an overall 
scale factor given by Q, the net charge on the island. 

We can now introduce a simpler charge non-conserving 
model. The majority of particles are neutral; further- 
more, the charged walkers tend to be repelled from the 
charged cluster, while the neutral, polarizable walkers 
are drawn to it. Hence one can expect that most of the 
growth will result from the aggregation of neutral parti- 
cles. Therefore, we n eglect the term due to the charged 
particles in Eq. ( |3.4[ ). The charge effects are taken into 
account simply by rescaling Q(t), As discussed above, the 
likely behavior of Q(t) is unknown, and dependent upon 
many competing effects. Here we will use the simplest 
possible rescaling rule: we will assume that Q{t) = A(t)a, 
where A(t) the overall area of the cluster and the average 
charge density a is assumed to be constant. 

We now have a sufficient set of ingredients for a com- 
plete growth model. That is, a given initial shape deter- 
mines the charge distribution and hence the electric field. 
The latter then allows a growth increment in time dt to 
be computed, yielding a new shape. With the new shape 
one then computes a new charge distribution and field, 
and so on. 



B. Analysis and Simulation 

Now we will show that, with the model as stated, a 
growing island will eventually transform from a compact 
to an elongated shape. Consider an island of elliptical 
shape, with principal axes a, b. If we compute the new 
boundary, using Eq. (3.4), the new shape will not be 
exactly elliptical; however, the ellipse's field may still be 
used as an approximation. It is then obvious that this 
new ellipse will have greater eccentricity than the original 
one if 



a + G(a,0)dt a 
b + G{0,b)dt > b 



(3.6) 



which (considering only growth from neutral walkers) is 
equivalent to 



G(a,0) 
G(M 



exp 



X 



E 2 (0) - E 2 (ir/2) /kT \ > a/b 



(3.7) 



The electrostatic field E{<f>) at distance 8 from the point 
(a cos 0,6 sin 0) is given by 



E 2 (<j>) 



28ab 



ill 2 i i 2-2 

(o cos <p + a sin 



'1/2 



(3.8) 



After introducing the average charge density a = Q/A, 
where A = nab is the total area of the cluster, Eq. (3.7) 
becomes 



exp 



-2 2 
X <7 IT 



2kTS 



(a-b) 



> 



(3.9) 



It is obvious from this equation that, if at some time 
during the growth the relationship ( |3.9p becomes valid, it 
will remain valid afterwards. Initially, the shape is com- 
pact (circular), i.e. a = b = R. However, due to the 
stochastic aspect of the growth there are always small 
variations in the radius. One can expect that the magni- 
tude of these variations is at least ~ 8, that is \a — b\ ~ 8. 
The transition from the compact to elongated shape oc- 
curs when the compact shape becomes unstable with re- 
spect to such variations. This happens when R reaches 
some critical value Rq, determined by 



1 



or 



8/Rq = ex P( % 2fcr ) 



2„2 



R w 28kT/ X a 2 TT 



(3.10) 



(3.11) 



Now, estimating a ~ ae/S^, \ ~ 2 x 10 24 cm 3 , 8 m ss 
8 ~ 10" 7 cm, a = 10%, and T w 300K, we get the 
estimate Rq ~ 0.2/im. Although this estimate is very 
rough, it is encouraging to see that Rq is smaller (by 
roughly an order of magnitude) than the size of seahorses 
observed in the experiments p4]Jl5[ |. 

Thus we find an instability of a compact, circular clus- 
ter to an elliptical form, when the compact cluster ex- 
ceeds a critical size. We have studied the growth of el- 
lipses numerically, and verified that the elongation decays 
for R < Rq; while for R > Rq the elongation persists and 
grows well beyond the linearized form, and in fact is am- 
plified (and 'pinched') by the resulting growth, to give 
two arms. We have also numerically tested instabilities 
to four arms. Here we find, for the parameters that we 
have explored, that the two-arm instability is dominant 
over four-arm instabilities. 

It is clear that, for most cases of interest, neither the 
charge density nor the electric field can be computed ana- 
lytically. Hence we need to implement our growth model 
numerically. This can be done as follows. The current 
boundary of the cluster is represented as a set of vertices 
{r v : v = . . . Ni> — 1}, connected by line segments. For 
each vertex we compute the potential energy of a walker 
near this point U(f v ); then a unit normal vector n v , and 
the displacement vector A„ = n v G(f v )A t are computed 
(where A t is a time-scaling constant). After all displace- 
ments are computed, all vertices are moved to their new 
positions. The only remaining problem is to compute the 
field distribution near the boundary. 

One can try to deal with this problem by solving 
straightforwardly Laplace's equation for the electrostatic 
potential. However, despite the fact that our problem 
is 2-dimensional, Laplace's equation for the electrostatic 
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potential is still 3-dimensional. Hence, to find the local 
electrostatic field, one would need to solve this 3D prob- 
lem for the full space. To take advantage of the lower 
dimension of the problem, we will use a surface charge 
method p9|-|32|], which we will outline below. 

Given the cluster C, we first compute the 2D charge 
density o~(x,y), and then find the electrostatic field. To 
find o~(x, j/), we introduce a grid {Di, i = 1, . . . , N g ; Di 6 
C}, forming a square lattice inside the cluster. Each 
square of this lattice is assigned a charge density <Xj (we 
will refer to these squares as elements). Now we replace 
the continuous equation for the electrostatic potential 
V(x, y) = constant for all (x, y) 6 C by the discrete equa- 
tion V(Di) = constant, which gives rise to the following 
set of linear algebraic equations: 



Uu = 2u;ln(5.44w//,) 



(3.16) 



Ng-l 

E 



(3.12) 



where i = 0, . . . , N g ~ 1, and the scaling constant on the 



right-hand side (here set to 1) can be chosen arbitrarily 
(that is, after the system is solved, we can rescale all 
the «7j). Uij is the electrostatic potential which would 
be induced at the point Di by the square element j, if 
the latter had a unit charge density. Therefore, one can 
approximate U^ as follows 



Aj 

U i:j = for i^j 



(3.13) 



\D* 



where Aj is the area of the j'-th element, and fj 
Di | . For most purposes the accuracy of this point-charge 
equation is sufficient, but one can also use more accurate 
expression for Uij which takes into account a quadrupole 
correction 



U, 



(1 



24^ 



The expression for the self-induced potential is 
U„ = 3.525 d 



(3.14) 



(3.15) 



where d is the element size (i.e. the grid step). 

Eq. ( 3.12| ) can be solved numerically and then the elec- 
tric field can be computed. However at this point we face 
a difficulty. As we have seen earlier, the field intensity 
E diverges near the edge of the cluster. In addition, the 
charge density itself diverges near the cluster's boundary. 
Thus we are trying to find small variations in a diverg- 
ing quantity, which itself depends on another diverging 
quantity. We can alleviate this problem by introducing 
line charges along the boundary. That is, in addition to 
charged square elements, we will assume that the bound- 
ary segments are also charged. This modification can 
be easily incorporated into Eq. ( 3.12| ); one needs only to 
distinguish between the self-potential Uu for square ele- 
ments, which is given by Eq. (3.15), and the self-potential 
Uu for line elements at the cluster's boundary, given by 



where w is the segment's width (i.e. border width) and 
li is the length of the segment. The accuracy can be 
further improved by using the exact field of a uniformly 
char ged r od rather than the point charge field given by 
Eq. fe.l3|) . 

As an even further simplification of the problem, we 
may use the linear (boundary) elements only. The ba- 
sis for neglecting the interior elements (which are exactly 
zero in 3D) is the fact that the charge density in 2D is the 
largest near the boundary, so that charges located close 
to the boundary give the main contribution to the electric 
field. This simplification not only reduces the number of 
elements to consider (and thus speeds the computation), 
but also makes it easy to compute the electric field in- 
tensity: to compute a field at distance 5 from element i, 
one can just use 



Ei = IhiOiwjb 



(3.17) 



The Boltzmann factor for a neutral particle in this field 
is then given by exp(\Ef / (kT)) = exp(7er^ 2 ), where 
7 = Axw 2 a 2 1 (S 2 kT) is the dimensionless interaction con- 
stant, and o~[ is t he d imensionless charge density, de- 
termined by Eq. ( |3.12 ) and by the rescaling condition 

£(oft)/EJi = i. 

We find this simplification to be very important, be- 
cause, when the charge density is modeled by a square 
grid (even if the grid is very fine), the nonuniformity 
in the charge density (arising from the finite grid step) 
causes noticeable (artifactual) fluctuations in electric 
field, which are then amplified by the growth of the is- 
land. In contrast (as will be seen below), if we locate 
all the charge on the boundary elements then the dis- 
cretization does not lead to such artifacts. Of course, this 
simplification results in a less realistic field, but we be- 
lieve that the main physical features of the system — the 
charge redistribution within the cluster, and the resulting 
branch repulsion and competition — are preserved. 

We add one further ingredient to our model, namely, 
the property (i): branching. The branching instability in 
growing forms is well known |}3] |55| ; it arises in diffusive 
problems as a result of competition between surface (in- 
coming) diffusion — which tends to favor new branches — 
and edge diffusion (along the boundary) — which tends 
to keep the cluster smooth. Such competition may be 
expressed in terms of a length scale Lq; this is essen- 
tially the cluster size, or branch length, at which edge 
diffusion can no longer maintain the smoothness of the 
growing form, and new branches form. Our continuum 
model does not explicitly represent the underlying diffu- 
sive processes. Hence we include a branching instability 
in a simple way as follows. When a branch reaches a size 
Lq (or after a certain number of growth steps — which is 
approximately the same criterion, for the fastest-growing 
branches) we force it to split exactly at the point of the 
fastest growth (as determined by the charge distribution) 
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into 2 branches. This branch splitting is just the simplest 
implementation of the fingering/tip splitting instability 
which is present in many growth phenomena || . This in- 
stability arises near the tips of fast advancing branches; 
and the faster the growth, the more likely the split to 
occur (3). For most diffusive problems one may expect 
branches to form also away from the fastest-growing tips. 
We have explored this possibility numerically, and find in 
fact that, if we introduce such branches away from the 
tips, they do not grow — due to the severe competition 
present in our model. Hence, in the following, we only 
describe cases where noise is introduced near the tips of 
the fastest-growing branches. 

Obviously, the model described thus far lacks any 
source of noise (apart from some tiny numerical noise, 
which is always inevitably present). This means (among 
other things) that all symmetries present in the initial 
shape of the cluster will be preserved during the growth. 
To explore the possibility of chiral symmetry breaking, 
then, we need to add some (very small) noise into the 
system — noise which breaks 2D inversion symmetry. 

Based on the considerations above, our typical starting 
condition for numerical studies is an ellipse, with some 
very small defects added to represent noise. As men- 
tioned above, the cases of interest, in which the defects 
induce further branching, occur when the noise is added 
near the tips — ie, near the long axis of the ellipse. As 
one test of the sensitivity of the model to chiral noise, 
we have used as a starting form an ellipse with 2 tiny 
defects, which are small dents placed near (but not at) 
the tips of the ellipse, along the long axis. The right 
dent is placed slightly above the x axis, and the left dent 
slightly below. Thus these defects break 2D inversion 
symmetry. Fig. [| shows the results after several genera- 
tions of branching. We note that the original defects are 
so tiny that they are almost invisible in Fig. [| however 
the original tiny chirality has been increased hugely by 
the resulting growth. We also see that Fig. |^ bears an 
obvious resemblance the experimentally observed "sea- 
horses". In fact, it has essentially the same geometric 
properties as a typical experimental seahorse: 1) only 2 
main arms; 2) these main arms are curved; 3) the cur- 
vature is correlated (that is, the two branches curve in 
opposite directions) ; 4) the outer edge of each main arm 
is covered by "fins" (small dead branches). 

It is important to test that such dramatic effects are 
actually implied by the dynamics of the model, rather 
than artifactual (ie, resulting from some defects in the 
implementation of the model). In fact, it is possible 
to "discover" a chiral symmetry breaking in the present 
model arising simply from the fact that a certain order- 
ing of points along the boundary (clockwise or counter- 
clockwise) is present in the numerical algorithm. Such an 
ordering can break the left /right symmetry of the numer- 
ical growth results unless sufficient care is exercised. We 
have discovered and removed such artifacts in our own 
algorithm. A good test is then a numerical simulation 
in which the starting shape is an ellipse with 2 defects 



which do not break the symmetry (ie, the defects are 
located precisely at the tips of the ellipse). The result 
of such a test is presented in Fig. ||: the original sym- 
metry is retained during subsequent growth. Hence we 
are confident that our numerical studies are free of any 
artifactual chiral symmetry breaking. 

Now let us consider in more detail how the observed 
curvature arises. Fig. [5] gives a more detailed picture of 
the growth and branching processes at one end of the 
island. After the first split (Fig. ||a), the new branches 
compete for further growth (as determined by the elec- 
tric field), and one of the new branches "dies" [ie, prac- 
tically stops growing (Fig. |^b)]. The winning branch 
eventually reaches the critical length and then splits 
again (Fig. |^c). Let us consider what happens when this 
branch splits in some detail. First, we note that the new 
branches "feel" one another. The electrostatic repulsion 
will cause the charges to redistribute, so that the points of 
fastest growth will no longer coincide with tips of the new 
branches. Instead, the fastest growth will occur at points 
of maximal local charge density (hence maximal |-B| 2 ), 
which arc farther away from each other than tips of the 
branches. This will force the branches to curve away from 
one another. If the split is originally symmetric, then the 
result of the competition between the new branches will 
be determined by the influence of all the other (already 
dead) branches. In particular, the main body (being 
closer to the lower branch) favors the upper new branch 
(branch 3), while the previous branch (branch 2) favors 
the lower new branch (branch 4). Branch 2 is smaller 
than the main body, but the main body is further away. 
Hence, if the parameter values are right (in particular, 
if the distance between splitting points is small enough, 
and dead branches are large enough), then the nearby 
dead branch has a dominant influence, and so branch 4 
will win the competition. Then, due to the exponential 
difference in the growth rate, branch 2 will stop growing. 
If this sequence is then repeated, the main arm curves 
further to the right at each branching. 

A similar process happens on the opposite end of the 
cluster. Hence we see how the present model leads to two 
main arms, each of which may be rather strongly curved. 
An 'S' shape (rather than a 'C or a '3' shape) results 
if the direction of curvature is the same for both main 
arms; that is, the two arms must show a correlation in 
curvature. This implies that one end of the growing body 
is significantly affected by the form of the distant end. 
We believe that this point is important: curvature itself, 
such as shown in Fig. |5[ may be ascribed to relatively local 
effects. However the correlation giving rise to consistent 
'S' shapes appears to require some form of long-ranged 
communication between the parts of the growing island. 
It is interesting to note here that, in a small minority of 
the experimentally observed islands, there are two arms 
whose curvature is in opposing directions — giving a '3' 
form rather than an 'S' [jig] . 

This correlation is not tested by our Fig. 0, since it is 
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built into the starting condition. We demonstrate that 
such a correlation does occur in the present model in 
Fig. 0. Here the starting condition is an ellipse with one 
(right) defect off axis (to break the symmetry), but the 
other one precisely on axis. We see that the growth dy- 
namics not only causes the right arm of the cluster to 
curve downwards, but also determines the outcome of 
the competition at the other end of the island. It is clear 
from the figure that the growth at this end is initially 
symmetric; yet subsequently, the upper arm dominates, 
purely due to influence from the other end of the cluster. 
We see that 'clockwise' growth at one end has sufficient 
influence to force clockwise growth at the other. Thus 
we find that dominance and curvature are correlated be- 
tween the two main arms, leading to strong chirality of 
the overall cluster. 

We have also examined the case that both defects are 
on the same side of the main axis. Here too we find that 
the two branches, while initially curving as expected from 
the local bias, after some growth 'feel' one another and so 
are repelled. We believe however that this starting con- 
dition is less realistic than those in Figs. ||and^, because 
the noise we introduce into our continuum model does not 
represent truly microscopic noise, but rather that noise 
that may be expected to grow beyond the microscopic 
scale. This is because the continuum model itself is not 
a truly microscopic model. Below (section [V) we will 



give evidence from our microscopic Monte Carlo studies 
that this continuum model does indeed capture impor- 
tant features of the large-scale behavior of a microscopic 
model with noise — in particular, the two-arm instability, 
and the repulsion between these two arms. 

It is entirely plausible that a "non-local" effect such 
as the correlation of curvature shown in Fig. ^| is more 
fragile than the more local effects shown in Fig. || This 
is true experimentally [jig] : as noted above, while the 
large majority of individual experimental clusters form 
an 'S' shape, some do not. Yet it is precisely this kind 
of long-distance correlation that gives rise to the broken 
inversion symmetry: in the absence of such correlation, 
half of the two-armed islands would have a '3' or 'C 
shape — and these shapes are not 2D chiral. 

Finally, we wish to demonstrate in yet another way 
that the dramatic symmetry breaking shown in Fig. |^ is 
a genuine outcome of the growth dynamics of our model: 
that is, the 'S' shape is not a necessary outcome of the 
growth model, yet neither does it require fine tuning of 
the model parameters. Starting with the same initial 
form as that leading to Fig. ^|, we have performed simu- 
lations for different values of the two model parameters: 
Lq (the branching length), and the electrostatic interac- 
tion constant 7. The results are presented in Fig. [?]. 

A number of conclusions may be drawn from Fig. |t]. 
We see clearly that branching is important for obtain- 
ing S shapes: in the bottom row of the figure, branch- 
ing is essentially absent, and the resulting forms amplify 
the broken symmetry of the initial condition only rather 
weakly. We also see that charge effects are equally impor- 



tant. They are weak in the left-hand column of Fig. [7]; 
and the resulting forms are uninteresting, regardless of 
the branching length. Finally, we note that, while both of 
these effects are important — in that they must be present 
in order that 'interesting' growth forms result — fine tun- 
ing of the two model parameters is not needed: we see 
clear two-armed 'S' shapes over a central region of the 
parameter space. 

In short: Fig. [?] — along with the previous figures — 
demonstrates rather clearly that growth dominated by 
(electrostatic effects + branching) can give rise to sponta- 
neous 'S' shapes which strongly break 2D inversion sym- 
metry, and that these forms are not simply modeling ar- 
tifacts arising either from fine tuning or from being built 
into the model by hand. 

We thus assert that our model possesses a genuine chi- 
ral growth instability; that is, it can drastically amplify 
any tiny perturbation that breaks the left/right sym- 
metry. Such 'drastic amplification' is of course what is 
meant by the term 'spontaneous symmetry breaking' — 
which is also an appropriate description. Hence the kind 
of pattern-forming behavior studied here falls into the 
same, quite rare class of behavior as that seen in phyl- 
lotaxis @J20|. 



IV. MONTE CARLO STUDIES 

We have performed a number of MC studies of aggre- 
gation under conditions where the aggregate is charged, 
and there are both charged and neutral walkers. The 
most obvious obstacle to overcome is the size of the clus- 
ters; as mentioned earlier, there are 10 6 -10 7 particles m a 
typical experimental cluster. Thus, in order to perform a 
successful MC simulation of these experiments, one must 
find a regime where the desired behavior (formation of 
S-like clusters) is achieved with smaller numbers of par- 
ticles. 

We have explored a number of MC models, involv- 
ing various combinations among the following choices: 
charge-conserving or not (as above); following each 
walker individually, or making (again) a quasiequilib- 
rium assumption for the spatial distribution of walkers; 
and treating the cluster as insulating or conducting. Our 
most encouraging results were obtained using the latter 
choice in each of the above three cases: that is, modeling 
only the neutral walkers, treated as a quasiequilibrium 
gas bonding to a conducting cluster. The resulting model 
is very like our continuum model in many ways — with the 
important difference that branching is now a part of the 
model dynamics, rather than being enforced as a premise 
of the model. 

With the above assumptions, the growth rule is as fol- 
lows. We perform the growth simulation with the under- 
lying 2D space discretized as a triangular lattice. Given a 
shape of the cluster (ie a set of occupied sites — our start- 
ing cluster is typically a single site), we first find all non- 
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occupied lattice sites adjacent to the cluster. Then, for 
all sites simultaneously, we compute the non-normalized 
probability of a new particle joining the cluster at this 
site according to 



P s = A exp[( X E 2 s /2 + K s E b )/kT] 



(4.1) 



where subscript s numbers sites along the cluster bound- 
ary, E s is the electric field strength at site s, A is a 
normalization constant (effectively, the time scaling con- 
stant), K s is the number of sites neighboring site s which 
are already occupied, and Eh is the bond energy. Thus, 
again we have a two-parameter model. The strength of 
electrostatic effects is set by the dimensionless \- The 
competition of surface and edge diffusion is contained in 
the bond energy E\,\ large Eh tends to keep the clus- 
ter smooth, while small E b allows small irregularities to 
grow. 

Using this MC model, we were able to reproduce some 
of the features of the seahorse clusters. Specifically, 
we found a clear two-arm instability for a significant 
area of the two-dimensional parameter space. We also 
found some tendency for the arms to curve, with small 
dead branches on the outer edges of the curves (Fig. ||). 
However — recognizing now that the particle number N 
represents in a sense a third dimension to be explored 
(with our present range of N < 2000) — we did not find a 
volume of this three-dimensional space in which the cur- 
vature appeared consistently, ft is however interesting 
to note that, where curvature did appear, there was an 
apparent tendency for the curvature of the two arms to 
be correlated. 

As mentioned above, in tests of our continuum model 
we found that there was a four-arm instability, which 
however was weaker than the two-arm instability over a 
range of growth parameters. One can also see this from 
Fig. for some parameters, one finds four dominant 
arms, and for others, two. Here we see the same com- 
petition in Fig. |§|: two arms have grown to dominance, 
yet there are two others which have not entirely stopped 
growing. We have also plotted out a two-dimensional 
'map' of MC growth patterns, depending on the elec- 
trostatic parameter \ an d the smoothness parameter 
Eh I kT. We find that this 'map' is qualitatively like that 
shown in Fig. [?]: compact shapes for small x an d large 
E b I kT (roughly equivalent to small 7 and large Lq in the 
continuum model); and two, relatively straight, arms for 
large x an d small Eb/kT. There are also some signifi- 
cant differences between the two maps. One difference 
is (again) that we find no region [in (x,Eb) space, at 
N ~ 2000] of consistent curvature of the two arms, giv- 
ing S shapes. Another is some tendency for three or six 
arms to appear; this is almost certainly an artifact of the 
underlying triangular lattice. 

Clearly these results support some important aspects 
of the picture gained from our other, less microscopic 
models, even as they fail to give 'seahorse' forms. The 
two-arm instability, with the tendency for the two arms 



to lie ps 180° apart, is clear; also the microscopic param- 
eters x an d Eb/kT give qualitatively the same growth 
behavior as do 7 and Lo in the continuum model. How- 
ever we have not found consistent S shapes from the MC 
models that we have explored. We believe that MC sim- 
ulations with larger TV are needed to test these ideas 
and results further, ft is clear that the kind of repeated 
branching and death shown in Fig. ^| requires a minimal 
number of microscopic particles in an aggregate before 
it can appear. Each dead branch which plays an impor- 
tant role in driving the curvature must be formed of some 
minimum number of particles; and there must be several 
generations of dead branches — and of course concomitant 
growth of the entire aggregate — for the curvature to be- 
come significant. We believe that our current results are 
below that threshold in N (assuming that it exists). Our 
current MC results also show too much dominance from 
both noise and the underlying lattice, ft is clear from 
Fig. H that artifacts from the underlying lattice are not 
negligible at the small scale of the figure. Also, micro- 
scopic noise is still large at this scale, compared to the 
other physical effects influencing the growth. Both of 
these effects will become less important at larger TV. 



V. GROWTH IN AN EXTERNAL 
ELECTROSTATIC FIELD 

The results that we have presented above represent an 
unusual form of symmetry breaking in the modeling of 
growth phenomena. While these studies were inspired 
by the same unusual symmetry breaking seen in the ICB 
experiments |]l4j-|l6|| , it is by no means certain that the 
mechanisms explored here are in fact responsible for the 
observed growth patterns. Our ideas have however mo- 
tivated further growth experiments, to be reported in 
detail in a separate publication |ll|. Here we will briefly 
describe the idea of the experiments, and the correspond- 
ing simulation studies that we have done, using our con- 
tinuum model. 

ICB growth experiments have been performed as in 
Refs. H] and [fL5[ , with the single change that there is 
imposed an external electrostatic field in the plane of the 
films. The motivation is to test whether electrostatic ef- 
fects are indeed important for the growth forms; and the 
results Q say that they are. Qualitatively, the in-plane 
field tends to give three effects: (a) the curvature of the 
main arms is reduced by the field, and even eliminated 
in a sufficiently strong field; (b) growth appears to be 
predominantly at one of the two main arms; (c) there 
is a weak tendency (which cannot be distinguished with 
certainty from zero) for the straightened clusters to grow 
in alignment with the external field. 

It is straightforward to include an external electric field 
in our continuum simulations, by adding another term 
to the equations (3.12) for the charge densities. A typi- 
cal result from such a simulation is presented in Fig. [|. 







The starting condition is the same as for Fig. | We 
see that the cluster has lost most of its curvature. Also 
there is some competition between the tendency to curve 
(which is favored by the starting defects) and the exter- 
nal field (which tends to align the growth along the field 
direction). The result of this competition is (roughly) a 
straight branch, growing at an angle to the field. These 
results show a good qualitative agreement with the exper- 
iments: (a) the curvature of the main arms is suppressed; 
(b) only one arm grows; and (c) there is some weak ten- 
dency for growth to follow the direction of external field, 
and at the same time, some reason to expect the growth 
to deviate from the orientation of the external field. 



VI. CONCLUSIONS 

Despite the obvious resemblance between the experi- 
mental results and our simulations, many questions need 
to be answered before one can claim that our model is 
directly related to the experiment. First, TCNQ crystals 
are not conducting (3^J|^]. However they are believed 
to have a large (about 8) dielectric constant Obvi- 
ously, the polarization charge distribution on the surface 
of a dielectric will be different from the charge distribu- 
tion on the surface of a conductor; but the larger the 
dielectric constant, the smaller this difference is. Sec- 
ond, our model relies crucially on the assumption of slow 
growth (this assumption leads to exponential differences 
in growth rates). It is not clear how close this assump- 
tion is to what really happens in the experiments. Third, 
there are very large uncertainties regarding the likely be- 
havior of Q(t), the charge on an island as a function of 
time. Finally, there is one feature in the experimental 
patterns that can almost certainly not be obtained from 
our electrostatic mechanism alone, namely, in the experi- 
mental seahorses the curvature is strong enough and pro- 
longed enough that the main arms commonly bend back 
to touch the central body of the island. We speculate 
that the addition of crystalline anisotropy effects at grain 
boundaries can yield this kind of behavior. 

One might argue that some other kind of long-range 
force — for example that coming from elastic effects — 
might produce an S shape. Of course we cannot rule 
out this possibility. However the physical effects ex- 
pected from elastic forces are rather different from those 
explored here, arising from electrostatic effects. Elastic 
stresses are not concentrated at growing tips, but rather 
at 'valleys' of the solid's boundary. Also, in our pic- 
ture the electrostatic effects are important beyond the 
boundary of the aggregate. Elastic forces can also extend 
through the substrate, beyond the boundary of a grow- 
ing island, during epitaxial growth on a crystalline sub- 
strate [[39); however such effects seem likely to be much 
smaller for an amorphous substrate such as that used in 
the ICB deposition experiments. Finally, it is not clear 
to us how elastic forces could give the broken symme- 



try seen experimentally, nor the strong sensitivity of the 
growth patterns to an external electric field Jl6| . We do 
note however that some combination of elastic and elec- 
trostatic effects may be needed to give rise to the high 
degree of 'bending back' of the main arms (noted above) 
seen in the experimental seahorses. 

We comment upon our modeling results. It remains 
to be demonstrated that all of the features reported here 
(two arms, branching, consistent correlated "deaths" of 
branches leading to curvature, and correlation between 
the curvatures of the two main arms) can be seen to oc- 
cur together in a single model, with none of them being 
imposed on the model by construction. Our continuum 
model shows all of these features but the branching it- 
self (which is uncontroversial) ; and our MC studies have 
shown hints (or more) of all of these features. We believe 
that all of our results, considered together, suggest rather 
strongly that the principal physical idea explored in our 
various models — trapped electric charge on a growing ag- 
gregate, leading to strongly enhanced growth in regions 
of strong electric field — can indeed lead to all of these fea- 
tures. However the final demonstration of this remains a 
challenge to future work in theory and modeling. 

In summary: we have proposed a novel mechanism for 
spontaneous chiral symmetry breaking during 2D aggre- 
gation. This mechanism does not rely on the presence 
of microscopic chirality; instead, it results from the exis- 
tence of a long-ranged electrostatic interaction in the sys- 
tem, due to a trapped charge on the growing aggregate. 
This interaction leads to a strong competition and repul- 
sion between growing branches, and — as we have shown 
here — can give rise to a strongly but spontaneously bro- 
ken 2D inversion symmetry. We have explored several 
different approaches to the simulation of growth in the 
presence of such an electrostatic interaction: tree models, 
continuum models, and Monte-Carlo simulations. Espe- 
cially encouraging results were obtained using the con- 
tinuum model; simulated growth often led to clusters 
which possess the same geometrical properties as the ex- 
perimentally observed "seahorses" |l4],|l5). The results 
of our simulations in the presence of an external elec- 
tric field are also in good qualitative agreement with the 
experimental results. Further theoretical work, and fur- 
ther experiments, are needed in order to clarify and test 
the connection between the ideas presented here and the 
"seahorse" experiments. However, using physical ideas 
motivated by the growth experiments, we have found a 
new class of 'electrostatic growth' models which spon- 
taneously break 2D inversion symmetry during growth; 
and we remain convinced that such spontaneous symme- 
try breaking occurs in the "seahorse" experiments, and 
so demands a theoretical model which does the same. 
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FIG. 1. The first two steps of a charged tree's growth. At 
every step, the tree selects a configuration with the lowest 
electrostatic energy. For the first step one of the rejected 
configurations is shown above the transition arrow. 

FIG. 2. The growth of a conducting tree. l\ — 0.7, lo = 0.5 
(arbitrary length units); 6 — 12°. In this case, growth accord- 
ing to an energy-minimization rule leads to an 'S' shape. 

FIG. 3. Sixty time steps in the growth of a charged island. 
The starting condition is an ellipse with two tiny defects, each 
placed just counterclockwise of the ellipse's symmetry axis. 
The tiny symmetry breaking present in these defects is dra- 
matically enhanced by the subsequent growth. 

FIG. 4. Growth of an ellipse with 2 tiny symmetrically lo- 
cated defects. Such a shape preserves its symmetry during 
the growth. Hence the symmetry breaking seen in the other 
figures is not artifactual. 

FIG. 5. Early stages in the growth of the right half of Fig. 3, 
viewed in detail. (a),(b) The broken symmetry in the original 
(tiny) defect leads to dominance of the lower branch, (c) The 
growing branch 1 splits again, (d) The nearby "dead" branch 
2 inhibits branch 3, so that again the lower branch dominates. 
In this way the 'memory' of the handedness (right or left) of 
the original small defect is maintained and amplified by the 
subsequent growth. 
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FIG. 6. Growth from an initial shape in which one defect 
(on the right) is off the symmetry axis, but the other (on the 
left) is on the symmetry axis. The broken symmetry on the 
right-hand side still has a large effect on the growth of the 
left-hand side — in fact, it causes the same (clockwise) branch 
to dominate there, even though the left-side defect is sym- 
metrically placed. 

FIG. 7. Growth forms for different values of the branching 
length Lo (increasing downwards) and the interaction con- 
stant 7 (increasing to the right). There is a region of the 
figure (roughly, the center) in which growth gives 'S' forms 
resembling those seen experimentally. However if the elec- 
trostatic effects are too weak (left side) or branching is too 
infrequent (bottom) then other forms result, displaying little 
or no chiral symmetry breaking. 

FIG. 8. A 2300-particle cluster grown using a Monte Carlo 
algorithm on a triangular lattice. The growth rule is described 
in the text; growth parameters are \ — 1-0 and Et/kT — 2.5. 
Although this cluster looks promising, we have not obtained 
such results consistently. It is probably necessary to use larger 
particle numbers to obtain such forms consistently. 

FIG. 9. Growth in the presence of an external electrostatic 
field, oriented horizontally to the right. Starting conditions 
and growth rule are otherwise as in Fig. 3. 
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